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We use Schwinger boson mean field theory to study the antiferromagnetic spin- 1/2 Heisenberg 
model on an anisotropic triangular lattice in the presence of a uniform external magnetic field. We 
calculate the field dependence of the spin incommensurability in the ordered spin spiral phase, and 
compare the results to the recent experiments in CS2CUCI4 by Coldea et al. (Phys. Rev. Lett. 86, 
1335 (2001)). 



The ground states of two-dimensionaU-jD) Heisenberg 
models continue to be of great interest. Em In this paper, 
we study the antiferromagnetic (AF) Heisenberg model 
on an anisotropic triangular-lattice in the presence of an 
external magnetic field along the z-axis, 



He 



(1) 



In the above equation, the summation runs over all the 
lattice sites i and their neighboring sites (i + S). We 
consider AF nearest neighbor spin-spin couplings, rep- 
resented by J and J' as shown in Fig. 1, with both J 
and J' > 0. In the absence of the field, this model is 
equivalent to a class ^-models recently considered by 
a number of authors .eIqEtEI The model includes several 
well known limiting cases. At J = 0, it is equivalent to 
2D square lattice model, whose ground state is a two- 
sublattice Neel phase. At J' = 0, it becomes a set of 
decoupled spin chains. At J' = J , it is reduced to the 
isotropic triangular-lattice model, where the ground state 
is a three-sublattice antiferromagnet. Experimentally, 
this model may be relevant to the insulating phase of the 
layered molecular crystals, ^(BEDT-TTF) 2 A,LI and 9- 
(BEDT-TTF) 2 RbZn(SCN) 4 ,el Our interest on this model 
is largely motivated by recent experiments on CS2CUCI4. 
That system is a^nuasi 2D S frustrated Heisenberg 

antiferromagnet Coldea et al.E3 have used the neutron 
scattering to study the ground state and the dynamics 
of the system in high magnetic fields. Among the obser- 
vations, these authors found that the incommensurate 
wave vector changes as the magnetic field increases, and 
the spiral spin density wave evolves into a fully saturated 
state. 

In this paper we apply the Schwinger boson mean field 
theory (MFT) to study the effect of magnetic field in 
the frustrated Heisenberg models. This method enables 
us to study incommensurate magnetic ordering in the 
quantum spin systems. The magnetic ordering is identi- 
fied as the Bose condensation of the Schwinger bosons, 
and the incommensuration of the ordering is determined 
by the wave vector of the condensed Schwinger bosons. 
In the absence of the magnetic field, the MFT pre- 




FIG. 1: The anisotropic triangle lattice: the bond along the 
b axis is J and the bond along the a axis is J'. 



diets three possible ground state: a two-sublattice Neel 
phase, a spiral spin state, and a spin liquid phase, simi- 
lar to the results obtained in the high temperature series 
expansions B In the presence of the field, we calculate the 
field dependence of the incommensuration in the spiral 
phase, and compare the results with the experimental ob- 
servation in CS2CUCI4 with good qualitative agreements. 

In terms of Schwinger bosons an and a^, the spin 
operators are expressed as 



a] 2 a i2 ) (2) 



with a local constraint at every site i given by a^an + 
a i2 a i2 = 1. As a standard method, we introduce a La- 
grangian multiplier field A.; to describe the constraint. 
The Hamiltonian of the system then becomes 



^ Matjftu + a i2 a i2 - 1) + 7T ^ Js 



(3) 



i,S 



where h = fj,B; fi, v = 1, 2, j = i + 5, and A. 
Oi/jO^v — Uivajfi is the spin singlet operator of bond (ij) 
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We note that operator Aij^ v is antisymmetric with re- 
spect to either the position (ij) or the indices (^v). On a 
bipartite lattice, a spin rotation by 7r on one of the sub- 
lattices transforms the spin-singlet bond operators into 
a symmetric operator with respect to the bond indices 
(ij). A mean field theory based on that transforma- 
tion has been developed by Auerbach and Arovas.Eil The 
method was extended to study the frustrated lattices by 
many authors .113 (For an overview of Schwinger boson 
theory, see Ref. [2] and references therein.) We intro- 
duce two types of mean fields: A A1!V (<5) = ^ (Ajj + s^ u ) 
and A = (Xj) where (• • • ) represents the thermodynamic 
average. The mean field Hamiltonian may be solved us- 
ing the conventional bosonic Bogliubov transformation 
as well as the Green function's method. 

The mean field Hamiltonian in Eq. (3) is diagonalized 

as 



Hmf — 



E 



(4) 



where the single boson spectra u> ^(k) = ui(k) ±h/2 with 



io{k) = ^A 2 -| 7 (£0| z and j(k) = £ 4 J 5 A 12 (6) sin(fc • 
8). (fi = ±). akf_i is a boson annihilation operator 
related to the original Schwinger bosons. £q/N\ = 
Etf, M „ ^(|A^( ( 5)| 2 + l/4)/2-2A. The free energy is given 
by (P = 1/ksT, with T the temperature) 



F = \ E Ml - ex P [-/3^(fc)] + ~ + £b- 

^ 1, L 7- L 

(5) 



The mean field Hamiltonian is solved together with the 
self-consistent equations for the two types of mean fields, 
which are given by 



E/(f^^yNK( fe ))+V2] 

fj, — ± 

= 2-60; 



E 



' ik Km(k-8)[n B (Lj^k)) + l/2} 



_ 2J S A(S) 



Js 



(2tt) 2 w(Jfe) 

- &o7(fe*) sin(fe* • S)/X 



(8a) 



(8b) 



The general features of the mean field solutions in 2D 
lattices are qualitatively given as below. At any finite 
T, b = 0. At T = 0, if min[w(fc)] ^ and b = 0, the 
ground state is a spin liquid, whose gap depends on the 
minimal value of the spectra u>(k). If min[w(fc)] = and 
&o is finite, the Schwinger bosons are condensed to the 
lowest energy state, and the system possesses the mag- 
netic long-range order. The ordering wave vector Q is 
determined by u)(k*) = 0. From the spin-spin correla- 
tions, 

X * a (q) = - Hm_ <T t S?(t)S«„(0)) (a=x,y,z), 



we have 

Xxx(q) = Xyy(q) = Xzz(q) 

1 v-^ / n B [uj(k)] 



A; 



u)(k) 



n B [uj(k + q)] + ~ 
u>(k + q) 



[\ 2 -l(kh(k + q]]-- 



(9) 



for h = 0, which indicates the mean field theory does 
not break SU(2) symmetry. In the thermodynamic limit, 
(N\ — > +00), the correlation functions are convergent 
except for Q = 2k*, 



1 v ^ A 



n B (Mk)) + 1/2] - 2; (6a) 



1 7(fc) 



sin(£; • 8) [n B (u}^(k)) + 1/2] = 2A(S). 



(6b) 



In the MFT, the magnetic orderinwpay be identified as 
the Schwinger boson condensation.£3 Below we examine 
the possible Schwinger boson condensation at the wave 
vector k — ±fc*, corresponding to the lowest energy of 
(D M (±fc*) — ► 0. We introduce a non- negative quantity, 



2A 



N A LJ.(k*) 



[n fl (w_(+fc*))+n fl (a;_(-fc*))], (7) 



such that the points k = ±k* in the integrals in Eqs.(|6 
and (6b) are taken into account separately, 



Xaa(Q) 



(10) 



which indicates that there exists long-range correlations 
with Q = 2k* whence b ^ 0. 

On the triangular-lattice (see Fig. 1), each site has 
six neighbors: ±ai, ±a 2 , ±(ai + a 2 ). It is convenient to 
write the wave vector k = (k±, k 2 ), with ki and k 2 the 
components of the vector along the directions of aiand 
0,2, respectively. The lattice constant a is set to be 1. 

We now consider the mean field solutions at T — 0. Let 
me first discuss the solutions in the absence of the field. 
At J' I J — 1 (the isotropic triangular lattice), fcp = 7r/3 
and u>(k*) = 0, 60 is finite indicating a magnetic long 
range order with the ordering wave vector Q(27r/3, 2tt/3). 
At J — (the square lattice), k* = tt/2 and oj(k*) = 0, 
and 60 > implying a Neel ordering at Q(tt,tt). The 
MFT in these two limiting cases is consistent with the 
known results. At J' = (decoupled one-dimensional 
chains), the MFT gives fep = 7r /4, (2k* along the chains), 
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FIG. 2: The critical field in unit J versus the ratio J'/ J. 



u)(k*) > 0, and b = 0, suggesting a spin gap state. The 
ID model is exactly, soluable, and the ground state is a 
gapless spin liquidE-3, although the static sin-spin corre- 
lation becomes strongest at Qf, = 7r.E£l. The discrepancy 
between the MFT and the exact solutions is primarily 
due to the neglect of the topological term in the MFT. 
For the general values of J'/ J, the MFT predicts three 
phases at ft, = 0: (1). a spin liquid phase at J'/J < 0.136; 

(2) . a spin spiral state at 0.136 < J'/J < 1.70, with 
the ordering wave vector bewteen (7r/2,7r/2) and (7r,7r); 

(3) . an antiferromagnet with the ordering wave vector 
Q = (n,ir) at J'/J > 1.70. We note that the phase di- 
agram of the model was studied previously by using a 
series expansion technique, liaaar-apin density wave and 
SP(N) mean field theory etcMffllid In the method of se- 
ries expansion, they found a Neel state persisting down 
to J'/J > 1.43, and predicts a spiral phase at small ra- 
tio of J'/J. These are consistent with our MFT. These 
authors also found a dimer phase between the Neel and 
spiral states. The dimer phase breaks translational in- 
variance and is not included in the present MFT. 

We now consider the soltuions at h ^ 0. In this 
case, the SU(2) rotational symmetry is broken, and 
uj-(k) < uj+(k). The Bose condensation criterion is given 
by uj-(k*) = u(k*) - h/2 = 0. We use the MFT to study 
the field dependent incommensurability in the spin spi- 
ral phase. The results are shown in Fig. 2, where the 
vector k* are plotted as a function of J'/ J at several val- 
ues of h. The main feature is as follows. 1) J'/J = 1 
is a stable fixed point, around which the external mag- 
netic field does not change the ordering wave vector; 2) 
At 0.136 < J'/J < 1, the ordering wave vector increases 
slightly as the field increases; 3) At J'/J > 1, the or- 
dering wave vector decreases as the field increases; 4) At 
J'/J < 0.136, the spin liquid may evolve into a spiral 
state, and becomes fully saturated as the field further 
increases. 

In the presence of the magnetic field, the spin z- 
component < S z >^ 0, and the ground state breaks 
the SU(2) invariance. At T = 0, the expectation value of 



FIG. 3: The ratio versus the wave vector k* at different 
fields are ploted. Inset: The wave vector near J'/J = 0.33 
at ten different fields are ploted. From left to right, h* — 
0.0, 0.1, - - ,1.0. 



S z is given by 

{Sf)= l -h*bo{k*,h*), (11) 

where h* = h/ (2A) is the dimensionless field. At h ^ 
0, (Sf) is finite, which indicates a polarized component 
along the field orientation z-axis. The static transverse 
susceptibilities at Q — 2k* are given by 

Xac(Q) = b&k*,h*) r _ 2 i 
N A 8 I 1 ' \ [ ' 

for a = x and y. Xaa{Q)/N A decreases as the field in- 
creases, and approaches to zero at h* = 1. b (k*,h*) 
can be calculated within the MFT. A special case is 
b (k*,h* — 1) = 1 corresponding to (Sf) = 1/2, or the 
spin full polarization. We have calculated the critical 
field [iB c (defined as the lowest field to induce the full 
spin polarization) as a function of J'/J. The results are 
plotted in Fig. 3. As we can see, /j,B c increases as J' 
increases for the fixed J. As J' increases, the spin cou- 
plings are strengthened, and it requires a higher field to 
polarize the spin. j— . 

Very recently, Coldea et al.Ej have reported the 
neutron scattering experiments on the antiferromagnet 
CS2CUCI4 in the high magnetic field. That system is 
a quasi-2D spin-1/2 quantum system in a triangular- 
lattice as shown in Fig.l. The spin-spin couplings are 
anisotropic with J'/J w 0.33. In the absence of the ex- 
ternal magnetic field, the spins are incommensurately or- 
dered and are aligned within the plane of the triangular- 
lattice. The latter may indicate a weak deviation from 
the Heisenberg model. Coldea et al. have studied the 
low temperature states of the system in the presence of 
in-plane as well as perperndicular magnetic fields. In the 
presence of the perpendicular fields, the states are found 
magnetically ordered with the varying incommensuration 
below a critical field, above which the system becomes 
fully spin polarized ferromagnet. In the presence of the 
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FIG. 4: The incommensuration e versus the external field at 
J'/J = 0.33. 



in-plane field, they have observed additional spin liquid 
phase between the incommensurate states and the ferro- 
magnetic phase. There have been theoretical efforts to 
understand their experimental results.cJ In the present 
paper, we have only considered the Heisenberg model in 
a uniform magnetic field. The predicted spin structure 
breaks SU(2) symmetry and shows the spin polarization 
along the field-direction. Such a spin structure is com- 
patible with the experiments in the perpendicular field, 
but incompatible with the in-plane field. Therefore, our 
MFT may be of relevance to the perpendicular field case 
in their experiments. 

To compare with the experiments, we define a quan- 
tity to describe the incommensuration, which is propor- 
tional to the wave vector deviation from the Neel state, 
e = 1/2 — Qi/tt, with Q\ the component of the spiral 
wave vector .aLong a\ direction. In Table I, we list the 
experimentally and theoretical values of e at field h* = 
and at the critical field h* = 1 for the two values of J'/ J. 
Also listed are the values of e at h* = calculated from 



the series expansion methodta. The agreements between 
the present MFT and the series expansion at h* — are 
very good. In Fig. 4, we plot the incommensuration rel- 
ative to the Neel state as a function of the external field 
for J'/J — 0.33. Our mean field results are in qualita- 
tive agreement with the experiments: as field increases, e 
also increases in the parameter space of interest. Quan- 
titatively, the theory predicts a weaker variation in the 
incommensuration than in the experiments. This dis- 
crepancy could be partly due to the neglect of the devi- 
ation of the physical system from the Heisenberg model 
in the theory. 



h* 


Exp.t£ 


J'/J = 0.404 


J'/J = 0.33 


SE0 


0.0 (e ) 


0.03 


0.047 


0.035 


0.028 


1.0 (e c ) 


0.053 


0.053 


0.0422 





Table I: The incommensuration to the Neel state along 
the b axis is listed, h* = means the absence the mag- 
netic field, and h* = 1 means the state is saturated fully. 
SE means the series expansion method- and the data is 
estimated from the work by Zheng et al.lj, which was also 
predicted by Manuel and CeccattoJa 

In summary, we have used a Schwinger boson mean 
field theory to study the antiferromagnetic Heisenberg 
model on an anisotropic triangular lattice in the presence 
of an external magnetic field. We calculate the magnetic 
field dependence of the incommensurability of the spin 
spiral phase. The theoretical results are compared well 
with the recent neutron scattering experiments. 

We would like to thank stimulating discussions with 
D. A. Tennant and R. Coldea on their experiments. This 
work was supported by a RGC grant of Hong Kong and a 
CRCG grant of the University of Hong Kong, and by the 
US DOE grant No. FG03-01ER45687, and by the Chi- 
nese Academy of Sciences. The authors also acknowledge 
ICTP at Trieste, Italy for its support and hospitality, 
where part of the work was initiated. 



1 S. Chakravarty, B. I. Halperin, and D. R. Nelson, Phys. 
Rev. B 39, 2344 (1989). 

2 A. Auerbach, Interacting Electrons and Quantum Mag- 
netism (Springer- Verlag, New York, 1994). 

3 Zheng Weihong, R. H. McKenzie, and R. R. P. Singh, Phys. 
Rev. B 59, 14367 (1999). 

4 J. Merino, R. H. McKenzie, J. B. Marston, and C. H. 
Chung, J. Phys.: Condens. Matter 11, 2965 (1999). 

5 L. O. Manuel and H. A. Ceccatto, Phys. Rev. B60, 9489 
(1999). 

6 C. H. Chung, J. B. Marston, and R. H. Mckenzie, J. Phys. 
Condens. Mater. 13, 5159 (2001). 

7 R. H. McKenzie, Comments Condens. Matter Phys. 18, 
309 (1998). 

8 J. Schmalian, Phys. Rev. Lett. 81, 4232 (1998); H. Kino 
and H. Kontani, J. Phys. Soc. Jpn. 67, 3691 (1998); H. 
Kino and T. Moriya, J. Phys. Soc. Jpn. 67, 3695 (1998); 
M. Votja and E. Dagotta, Phys. Rev. B59, 713 (1999). 



9 R. Coldea, D. A. Tennant, R. A. Cowley, D. F. McMorrow, 
B. Dorner, and Z. Tylczynski, Phys. Rev. Lett. 79, 151 
(1997). 

R. Coldea, D. A. Tennant, A. M. Tsvelik, and Z. Tylczyn- 
ski, Phys. Rev. Rev. Lett. 86, 1335 (2001). 

1 A. Auerbach and D. P. Arovas, Phys. Rev. Lett. 61 , 617 
(1988); D. P. Arovas and A. Auerbach, Phys. Rev. B 38, 
316(1988). 

2 I. Richey and P. Coleman, J. Phys. Condens. Mater. 2, 
9227 (1990); D. Yoshioka and J. Miyazaki, J. Phys. Soc. 
Jpn. 60, 614 (1991); A. Mattsson, Phys. Rev. B 51, 11574 
(1995). 

3 S. Sarker, C. Jayaprakash, H. R. Krishnamurthy, M. Ma, 
Phys. Rev. B 40, 5028 (1989). 

4 H. Bethe, Z. Phys. 71, 205 (1931). 

5 S. Q. Shen, Inter. J. Mod. Phys. 12, 709 (1998). 

6 M. Bocquet et al., Phys. Rev. B 64, 094425 (2001). 



